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ABSTRACT 


When  moving  air  which  is  in  equilibrium  with  an 
underlying  surface  of  a  given  constant  roughness  abruptly 
encounters  a  surface  with  a  different  roughness  length  a 
new  equilibrium  condition  is  set  up  gradually.  Several 
theories  have  been  proposed  to  explain  the  process  of 
transition  under  neutral  conditions.  The  author  s 
experimental  wind  profiles  are  presented  and  compared  to 
the  theoretical  profiles  obtained  from  a  modified  form  of 
Townsend’s  theory  of  self-preserving  flow  in  a  turbulent 
boundary  layer.  Discrepancies  between  theory  and  observation 
are  discussed  and  suggestions  are  presented  for  obtaining 


definitive  results. 
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INTRODUCTION 


The  study  of  the  effects  of  the  earth’s  surface  on 
the  atmosphere  is  of  prime  importance  for  future  advances  in 
numerical  weather  prediction.  Over  a  period  of  1  or  2  days 
energy  transfer  changes  between  the  surface  and  atmosphere  can 
be  neglected  but  for  longer  periods  such  changes  must  be 
included.  A  great  deal  of  work  has  been  done  on  vertical 
fluxes  of  heat,  momentum,  and  moisture,  and  on  the  wind  profile 
in  the  various  regimes  of  stability  particularly  over  uniform 
terrain.  In  this  paper,  I  shall  consider  the  effect  of  a  single 
change  of  surface  roughness  on  the  wind  in  the  surface  boundary 
layer. 

Suppose  wind  has  been  blowing  over  a  region  of  uniform 
roughness  length  for  a  long  period  of  time  and  that  the  atmosphere 
is  in  neutral  equilibrium.  Then  the  wind  and  turbulence  will  not 
be  changing  with  x,  the  downwind  distance  coordinate,  and  the 
wind  is  said  to  be  in  equilibrium  with  the  surface. 

Suppose  now  that  the  surface  abruptly  changes  to  one  with 
an  increased  roughness  length.  For  example,  the  wind  blows  from 
an  extensive  flat  sandy  area  to  a  region  of  grass.  The  first 
effect  on  the  wind  is  increased  friction  in  the  thin  layer 
immediately  adjacent  to  the  surface  and  thus  the  lowest  layer 


? 
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undergoes  a  deceleration.  The  velocity  gradient  in  the  vertical 
increases,  and  this  in  turn  increases  the  turbulence.  The 
increased  turbulence  feeds  more  momentum  into  the  lower  layers  to 
partially  overcome  the  momentum  lost.  Gradually  successively 
thicker  layers  of  the  surface  boundary  layer  are  affected,  and 
eventually  the  wind  again  will  be  in  equilibrium  with  the  under¬ 
lying  surface,  but  with  a  slightly  lower  mean  velocity  at  any 
given  level  and  with  an  increased  surface  shear  stress.  Also 
from  the  equation  of  continuity  the  overall  retardation  will 
induce  a  vertical  velocity. 

Suppose  we  measure  the  velocity  profile  under  conditions 
of  neutral  stability  just  ahead  of  the  surface  discontinuity  and 
at  several  places  downstream  as  in  Figure  1. 


Fig.  1 
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It  would  be  informative  to  see  what  qualitative  changes 

in  the  profile  would  occur.  It  is  convenient  to  graph  the 

profiles  with  log  z  against  u  rather  that  z  against  u. 

Let  us  assume  that  the  surface  roughnesses  are  z  and  z. 

o  1 

upstream  and  downstream,  respectively.  Because  of  the  assumed 
stability  let  us  assume  that  both  the  original  and  the  new 

profiles  are  described  by  the  logarithmic  profile  (at  least  in 

u*  u* 

2  *1  2 

the  lowest  layer),  e.g.,  u  =  ~r~  In  —  ,  x  <  0;  u  =  _  In  —  ,  x  >  0. 

K  Zo  k  Z1 

Then  if  z^  >  z^  the  following  set  of  profiles  might  be  expected. 


/  03  2 


U 


Z, 


The  original  profile  is  a)  and  lines  b) ,  c) ,  d  are  the 
profiles  modified  farther  and  farther  downstream.  The  point  of 
intersection  of  each  profile  with  the  original  profile  serves  to 
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mark  the  top  of  the  modified  layer  or  internal  boundary  layer  ,  h  .  It 
is  observed  from  experiments  that  the  depth  of  the  internal  boundary 
layer  indeed  grows  with  x.  If  z^  >  z  ,  the  following  profiles 


might  be  expected: 


,  *2  o  ,  , 

Friction  velocity  is  defined  as  u  -  where  t  is 

P  0 

surface  shear  stress  and  P  is  the  density.  In  this  paper  the 

*2 

kinematic  shear  stress  will  be  used.  It  is  defined  as  u  =  tq. 


The  friction  velocity  is  a  function  of  x.  On  the  graphs 

the  slope  m  =  .  In  other  words  the  friction  velocity  and 

u 

therefore  the  square  root  of  the  surface  stress  is  inversely 


5. 


proportional  to  the  slope.  The  constant  k  is  approximately 
equal  to  0.4  and  is  called  Von  Kerman* s  constant. 

Now  certain  qualitative  conclusions  can  be  drawn  and 
questions  asked: 

(1)  The  boundary  layer  grows  with  x.  What  is  the  functional 
relationship  between  h  and  x? 

(2)  (a)  In  the  case  of  a  smooth  surface  changing  to  a  rough 

surface  the  friction  velocity  is  at  a  maximum  close 
to  the  discontinuity  and  gradually  approaches  a  new 
equilibrium  value  greater  than  the  original  value  but 
much  smaller  than  the  maximum. 

(b)  In  the  case  of  a  rough-to-smooth  transition  the 

friction  velocity  is  at  a  minimum  immediately  behind 
the  discontinuity  and  then  gradually  approaches  a  new 
equilibrium  value  less  than  the  original  value.  How 
does  the  surface  stress  vary  with  x? 

(3)  This  simple  model  is  unrealistic  in  that  the  shear  stress  is 

assumed  constant  with  height  in  the  modified  layer.  It  must 
vary  smoothly  from  the  surface  value  to  the  value  at  the 

top  of  the  layer  which  is  x  .  What  will  be  the  shape  of  the 
velocity  profiles  allowing  for  the  smooth  vertical  variation 
of  stress?  Will  a  simple  linear  variation  of  stress  with 
height  account  for  the  shape  of  the  velocity  profile  or  is 

a  more  involved  relationship  required? 


* 

, 


The  following  theories  will  be  examined  in  Chapter  I: 
those  of  Elliott  (1958) ,  Taylor  (1962) ,  Panofsky  and  Townsend  (1965) , 
a  modification  to  Townsend's  work  by  Blom  and  Wartena  (1969)  and 
a  numerical  model  by  Smith  (1967).  Chapter  II  will  deal  with  my 
experimental  results  and  Chapter  III  will  offer  suggestions  for 
future  experiments. 
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CHAPTER  I 


One  of  the  first  theoretical  studies  of  the  effects  of 
a  change  in  roughness  was  made  by  William  P.  Elliott  in  1958. 

He  assumed,  as  have  subsequent  workers,  an  atmosphere  in 
neutral  equilibrium.  Then  he  assumed  that,  downwind  from  the 
change  in  roughness,  the  lowest  layer  was  in  equilibrium  with 
the  new  surface  up  to  a  height  h(x)  and  was  characterized  by 
the  logarithmic  wind  profile.  This  layer  was  called  the  internal 
boundary  layer.  Above  the  internal  boundary  layer  it  was  assumed 
that  the  air  was  as  yet  unaffected  by  the  change  in  roughness. 
Elliott  employed  Von  Reiman' s  integral  theorem  to  get  an 
expression  for  h.  He  assumed  logarithmic  wind  profiles  for  u 
in  the  internal  boundary  layer. 

The  chief  difficulty  encountered  with  this  approach, 
besides  the  clumsy  solution,  is  the  implied  discontinuity  in 
stress  at  the  boundary  surface  z  =  h. 

In  order  to  overcome  Elliott's  unrealistic  discontinuity 
of  stress  at  the  interface,  Panofsky  and  Townsend  used  a  log- 
linear  velocity  profile  below  the  interface 

* 

U  -  ^  E  (1-s)  In  +  3  £]  (1) 

k  z^  h 

*  * 

u  -u 

where  s  =  — ^ -  and  h  is  the  height  of  the  interface. 

u  * 
o 


j 
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This  can  be  rewritten: 


*  * 
where  uq  represents  friction  velocity  at  x  <  0,  and 

represents  the  friction  velocity  at  x  >  0.  Note  that  when 

z  <<  h;  the  distribution  approaches  the  logarithmic  one. 


Panofsky  and  Townsend’s  final  result  gives  h 
implicitly  as  a  function  of  x  in  the  following  equation: 


7  12 

4k3 (£-£,  )  =  n[ln  n  -  5  -  r  M  +  - ; 

°  Z  lnn-l-^M 


+ 


(In  n  -  1  - 


(2) 


+  terms  of  order 


[Inn  -1  -j  M]  3]. 


where  5  -  x/z  ,  n  =  h/z  ,  M  =  In  z  /z  .  The  two  profiles 
o  o  ox 

are  very  similar.  Differences  result  from  the  assumed  linear 
variation  of  stress  from  the  surface  to  the  top  of  the  boundary 
layer  which  leads  to : 

(a)  no  discontinuity  of  the  velocity  gradient 

(b)  no  discontinuity  of  the  shear  stress 

(c)  the  resulting  height  of  the  boundary  layer  is  greater 
(mainly  because  of  a  change  in  definition  of  the  location 
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of  the  upper  surface) . 

In  1965  Townsend  published  two  important  papers  dealing 
with  the  surface  boundary  layer  and  employing  the  ideas  of 
similarity  theory.  His  first  paper  (Townsend  1965a)  developed 
the  fundamental  concepts  of  self-preserving  flow  while  in  his 
second  paper  (1965b)  he  applied  these  ideas  to  particular 
special  cases  such  as  flow  over  a  hedge,  flow  over  a  fence  and 
flow  across  a  change  in  roughness.  He  examined  the  effects  on 
the  temperature  and  moisture  profiles  as  well  as  on  the  velocity 
profile.  It  was  shown  that  changes  in  the  flow  such  as  those 
caused  by  a  change  in  roughness  are  substantially  independent  of 
any  special  assumption  about  the  interaction  between  turbulence 
and  the  mean  flow. 

Self-preserving  flow  implies  that  lateral  distributions 
of  various  quantities  have  the  same  form  at  all  distances  from 
the  origin,  differing  only  in  common  scales  of  velocity  and  length. 
In  the  cases  considered  by  Townsend  it  was  assumed  that  the 
deviation  of  a  flow  quantity  from  its  value  in  the  undisturbed 
flow  was  self  preserving.  It  was  shown  in  his  original  paper  that 
self-preserving  development  in  the  atmospheric  boundary  layer 
occurred  if: 

(1)  The  velocity  defect  ratio  |  u-u.J /u^  was  small  except 

for  small  z/£  . 

s 

(2)  The  variation  of  the  roughness  with  x  was  small. 
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(3)  log  l  It.  was  large, 
s  i 


The  quantity 


will  be  defined  later. 


As  in  the  previous  studies  Townsend  started  with  a  wind 

blowing  from  an  infinite  plain  of  roughness  zq  at  x  <  0  onto 

a  field  with  roughness  length  z^  at  x  >  0.  The  atmosphere  was 

assumed  to  be  neutral  with  the  wind  profile  given  by  u  =  In  . 

o 

There  are  two  contributions  to  the  modification  of  the 
wind  profile.  One  is  due  to  the  acceleration  of  particles  of  air  V(z) 
and  the  other  is  the  vertical  displacement  of  the  streamlines  <$(z). 

The  profile  for  x  >  0  he  wrote  as: 


u  =  u  +  V(z)  -  —  6(z) 
o  kz 


(3) 


where  is  the  original  unmodified  wind  profile. 

Near  the  surface  the  profile  immediately  adapted  to  the  new 
friction  velocity.  Thus 


u  "  ¥  (Ti(x))i5ln  t 


served  as  a  lower  boundary  condition. 


To  find  a  solution  for  V(z)  Townsend  assumed  a  self¬ 
preserving  flow  of  the  form 


u 

V  =  -jjy-  f(n) 


(4) 
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where  u  =  u  (x)  is  a  velocity  scale,  n  =  z /£  ,  £  =  £  (x) 

s  s  s  s  s 

is  a  length  scale,  and  f(n)  is  a  universal  function. 

The  shear  stress  was  also  assumed  to  be  self-preserving 
and  of  the  form: 


T 


+  t  F(n) 
s 


(5) 


where  t  =  i.  -  u  and  F(n) 
s  1  o 

Taking  the  equation  of  motion 


is  another  universal  function. 


u 


3_u 

3x 


+  W 


a  u 

3  z 


ax 

dZ 


and  substituting  equation  (3)  and  (4)  he  obtained 


_  ,df.  dF 

■ n  W  ■ 


(6) 


To  obtain  an  explicit  form  for  F(n)  and  f(n)  it  was  necessary 
to  make  an  assumption  regarding  the  interaction  of  the  mean 
flow  and  turbulence.  Townsend  used  the  mixing-length  momentum 
transfer  relation 


du  _  x_ 
dz  kz 


which  led  to 
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df  =  F 
dn  n 


Combining  (6)  and  (7)  yielded  the  results: 


co  _x 

F(n)  c  e“n;  f(n)  =  ~  /  “ “  dx. 

n 


(7) 


(8) 


Townsend  determined  the  velocity  and  length  scales  ug  and 

£  using  the  lower  boundary  condition  together  with  the 
s 

universal  functions  f(ri)  and  F(ri)  .  He  arrived  finally  at 
the  result 


U  =  U - [/  -  dx(l+P  )  +  V  (1-e  n)n  ]  (9) 

o  k  J  x  1  1 

n 


where 


P,  =  [In  £  / z  -  C  ] 

1  so  o 

£  [In  £  /z  -  M  -  1]  =  2k2x 
s  si 

M  =  In  z  /z. 

o  i 

u  =  -  Mu  {[In  £  /z  -  C][l  +  P.]} 
s  O  si  1 

C  -  .577  C  =  1.577 
o 


For  the  surface  shear  stress  he  found 
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2  & 

i  *  =  u  +  u  (1+P.)  (10) 

1  o  s  1 

In  summary  Townsend  showed  that  flow  disturbances 

caused  by  a  change  in  roughness  are  self-preserving  if  certain 

conditions  are  met  in  the  atmosphere.  The  requirement  that  the 

distribution  of  velocity  should  be  logarithmic  near  the  surface 

makes  the  prediction  of  surface  shear  stress  nearly  independent 

of  the  exact  nature  of  the  turbulent  transfer  process  and  the 

velocity  profile  is  determined  within  narrow  limits  by  the  surface 

fluxes.  To  provide  an  explicit  profile,  some  assumption  about  the 

nature  of  the  interaction  between  the  mean  flow  and  turbulence 

is  required  in  order  to  derive  an  expression  for  f(n). 

°°  -x 

Townsend  used  f(q)  =  -  /  e  /x  dx  and  showed  that  Elliott 

n 

implicity  used  f(n)  =  In  n  for  q  <  1  and  f(q)  -  0  for  n  >  1 

while  the  "Panofsky  and  Townsend"  profile  used 

f(q)  =  In  jn  +  (1-  -j)  f°r  n  <  2  and  f(n)  =  0  for  q  >  2. 

In  a  recent  paper  Blom  and  Wartena  (1969)  made  a 

slight  modification  to  the  work  of  Townsend.  In  this  paper 

the  authors  noted  an  inconsistency  between  Townsend's  profile 

and  his  lower  boundary  condition.  This  inconsistency  arose 

because  of  the  first  order  approximation  Townsend  used  in 

evaluating  the  integral  /  V  dp .  He  used  the  approximate 

o 

solution  qV.  This  integral  arose  out  of  the  derivation  of  an 
expression  between  V(z)  and  6(z).  Blom  and  Wartena  were 
able  to  give  an  exact  expression  for  the  integral  using  the 


! 
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-x 


—  fl  6 

functions  F(q)  =  e  and  f(n)  =  /  -  dx  to  obtain 


x 


n 


n  u  n  00  -x  u  °°  —x 

'  Vdn  =  ~  /  f  — —  dx  =  -7^-  [r)  /  — —  dx  -  1  +  e  ]  .  (11) 


O 


k  "  x 

o  n 


k  "  "  x 

n 


For  small  n  this  gives 


nu. 


/  Vdn  =  nv  -  —  . 


(12) 


This  led  to  the  following  expression  for  the  velocity  scale: 


u  =  -  Mu* (1+Pn )  1 [ In  £  /z.  -  C  +  (1+P.1)  X]  1.  (13) 

s  o  1  s  1  1 


The  new  velocity  scale  now  allows  a  domain  to  exist  which 
satisfies  the  lower  boundary  condition. 

At  the  present  time  the  Townsend  approach,  as  modified 
by  Blom  and  Wartena,  is  the  most  physically  secure  theory 
available  which  treats  the  problem  of  change  in  roughness  in  the 
neutral  atmosphere. 

Another  interesting  solution  was  developed  by  Smith  (1967) . 
He  used  the  momentum  equation 


*2 

dn  _  _1  _3t  9u 
3x  P  3z  dz 
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*2  du 

the  momentum  flux  equation  (u  =  K  —  ,  where  K  is  the 

o  Z 

diffusivity) ;  and  an  empirical  equation  for  the  rate  of  change  of  K 


,  3K  .  u  3u 

(u  ^  =  M  ~  n 

u 


K 

3z 


A  is  an  empirical  constant  which  must  be  found  from  experimental 
data.  These  three  equations  were  solved  numerically  with  the 
lower  boundary  condition  K  =  0  at  the  ground.  The  resulting 
profiles  of  u  and  u  appear  realistic,  with  a  new  equilibrium 
regime  being  established,  given  sufficient  fetch  and  with  a  return 
to  the  original  profile  when  the  roughness  returns  to  its 
original  value. 

One  advantage  of  this  approach  is  that  the  new  roughness 

length  can  be  allowed  to  vary  periodically  about  a  mean  value. 

Smith  found  that  the  effective  roughness  length  depended 

essentially  on  the  period  of  oscillation  and,  to  a  smaller 

* 

extent,  upon  its  amplitude.  The  equilibrium  value  of  u  was 

* 

found  to  be  a  few  percent  smaller  than  the  corresponding  u 
for  a  constant  roughness. 

Taylor  in  1962  arrived  at  a  different  conclusion 
regarding  the  variation  of  surface  shear  stress  with  distance. 
Whereas  the  other  theories  predict  that  the  surface  stress 
gradually  approaches  an  equilibrium  value  far  downstream. 


! 
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Taylor  states  that  the  new  equilibrium  value  shoud  be  achieved 
almost  instantly.  Wind  tunnel  measurements  conducted  by  Taylor 
tended  to  support  this  idea. 

In  view  of  the  lack  of  agreement  concerning  the  downwind 
shear  stress  future  experimental  investigations  should  attempt 
to  determine  it  accurately. 


. 


CHAPTER  II 


THE  SITE 

In  order  to  obtain  as  large  a  change  in  surface  roughness 
as  possible  itT  was  decided  to  use  a  water-to-land  discontinuity. 
This  would  emphasize  the  effect  of  roughness  and  tend  to 
minimize  elevation,  stability,  and  possible  upwind  effects. 

The  requirements  for  a  suitable  location  were  as  follows: 

(1)  gentle  slope  away  from  the  water's  edge 

(2)  long  fetch  over  the  water 

(3)  flat  land  surrounding  the  lake  with  long  fetches  over  uniform 
grass 

(4)  easy  access. 

Only  one  lake  was  found  near  Suf field,  where  the  work 
was  carried  out,  which  would  satisfy  these  requirements.  It  was 
called  Scott's  Slough  and  was  located  two  miles  south  of  the 
research  station. 

The  land  surrounding  the  lake  was  very  flat  to  the  west 
and  north.  To  the  east  and  south  were  gentle  rolling  hills.  The 
local  relief  was  such  that  only  east,  south,  and  northwest  winds 
could  be  used.  The  reasons  the  other  directions  were  unsuitable 
were : 
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(1)  Insufficient  stretch  of  uniform  grass  in  many  downwind  areas, 

(2)  A  wide,  soft  mud  flat  on  the  south  shore, 

(3)  Weeds  growing  out  of  the  water  along  the  s'outhwest  shore, 

(4)  A  large  island  blocking  the  east  shore, 

(5)  Part  of  the  south  shore  rising  rapidly  away  from  the  water. 

DESCRIPTION  OF  THE  EXPERIMENTAL  AREA 
AREA  A 

Area  A  was  located  on  a  peninsula  and  was  used  when 
the  wind  was  roughly  from  300  degrees.  The  fetch  was  approximately 
280  meters.  The  surface  consisted  of  3  meters  of  mud,  followed 
by  16  meters  of  tall  (70  cm)  fox  tail.  After  the  fox  tail 
was  another  16  meters  of  shorter  (30  cm)  grass  before  the 
grass  became  the  short  blue  prairie  grass. 

AREAS  B  AND  C 

These  areas  were  used  with  winds  from  90  degrees. 

The  fetch  was  about  300  meters.  The  surface  consisted  of 
3  meters  of  mud  which  changed  abruptly  to  50-60  cm  foxtail 
that  extended  for  a  distance  of  30  meters.  The  next  40-45  meters 
consisted  of  shorter  (20  cm)  denser  grass.  Point  C  has  a  more 
uniform  density  and  height  of  grass  than  point  B. 
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f<3-  4  Scott's  Slough 
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AREAS  D  AND  E 

These  areas  were  used  with  180  to  230  degree  winds. 
The  fetch  was  350  meters.  The  surface  was  the  best,  being  mud 
for  4  meters  before  becoming  very  tall  uniform  foxtail  for 
60  meters  at  E  and  40  meters  at  D  with  shorter  grass  beyond. 
In  all  areas  the  transition  from  one  type  of  grass  to  another  took 
place  along  a  very  distinct  line.  The  mud  f lat-to-grass  transition 
was  also  very  abrupt. 


INSTRUMENTATION 

We  used  two  8-m  masts  and  one  16-m  portable  tower 
permanently  mounted  on  a  trailer.  The  masts  were  constructed  of 
2-in.  diameter  pipe  with  horizontal  supports  placed  such  that 
the  anemometer  cups  were  located  0.5,  1.0,  2.0,  4.0,  and  8.0  meters 
above  the  ground.  The  portable  tower  had  been  damaged  in  a  high 
wind  so  that  it  was  not  possible  to  use  it  to  obtain  the 
16.0  meter  wind  except  on  two  early  trials.  It  held  the  anemometers 
at  the  same  heights. 

The  anemometers  were  Shepherd-Cassella  cup  anemometers. 

Their  contacts  were  designed  to  close  momentarily  for  every 
1/720  mile  of  wind  which  passed.  The  anemometers  were 
connected  to  electric  counters  with  multi-conductor  cable. 


The  usual  procedure  was  to  place  one  mast  as  close  as 
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practicable  to  the  water  on  the  mud  flat.  The  mud  under  water 
was  very  soft  and  unsuitable  for  anchoring  a  guy  rope.  As  a 
consequence  the  mast  had  to  be  two  or  three  meters  from  the 
water's  edge.  However  there  was  probably  very  little  change  in 
roughness  between  the  water  and  the  mud  compared  with  the 
transition  from  the  mud  to  the  grass. 

The  second  mast  and  the  tower  were  positioned  downwind 
from  the  first  mast.  Wind  direction  was  determined  by  means  of 
a  portable  wind  vane  at  2  meters  above  the  ground.  It  was  set 
up  for  each  trial  with  a  compass.  Distances  between  masts  were 
measured  with  a  surveyor's  tape.  Air  and  water  temperatures  were 
measured  on  several  days  with  standard  mercury  thermometers. 

The  counters  and  vehicle  were  positioned  well  to  one 
side  of  the  line  of  towers,  (usually  more  than  25  meters).  Every 
effort  was  made  to  minimize  disturbance  of  the  grass  while 
setting  up  towers.  Estimates  of  the  fetches  were  obtained  from 
aerial  photographs. 

The  equipment  was  easily  portable  and  could  be 
transported  to  the  site  and  set  up  by  two  people  in  two  hours. 

The  procedure  followed  each  day  was  to  obtain  the  weather 
forecast  in  the  morning  from  the  Meteorology  Section.  If  the 
wind  was  forecast  to  be  from  a  suitable  direction  later  in  the 
day  the  towers  were  set  for  the  forecast  wind  direction  and  a 
test  was  carried  out  if  the  wind  cooperated.  This  procedure 
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yielded  useful  data  on  nine  days.  The  most  common 
reason  for  failure  to  obtain  more  data  was  the  lack  of  good 
pressure  gradients  during  much  of  the  summer  of  1968.  As  a 
consequence  the  winds  were  mostly  light  and  variable. 


THE  DATA 


The  -counters  which  recorded  wind  speed  were  stopped 
and  read  every  15  minutes.  The  total  amount  of  useful  data 
from  each  trial  averaged  about  2  hours.  The  results  were 
divided  into  two  groups.  One  group  consisted  of  data  obtained  on 
sunny  days  (assumed  to  be  unstable) .  The  other  group  was 
obtained  under  overcast  conditions  (assumed  to  be  near  neutral) . 

On  5  of  these  occasions  the  towers  were  placed  at  3,  33,  and 
63  to  78  meters.  An  overall  mean  profile  was  plotted  for 
these  cases  (Figure  5)  .  They  consisted  of  3  unstable  cases 
and  2  neutral  cases.  Mean  profiles  were  also  plotted  for  the 
neutral  cases  (Figure  £>)  and  the  unstable  cases  (Figure  7), 
separately. 

Two  trials  were  conducted  with  towers  located  at  3,  13, 
and  23  meters  from  the  edge  of  the  water.  One  was  neutral  and 
the  other  was  unstable  (Figure  S) . 

From  these  graphs  it  is  possible  to  measure  the  roughness 
lengths  and  the  friction  velocities.  Knowing  these  parameters  it 
is  possible  to  make  some  comparisons  with  theory. 
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The  profiles  obtained  from  each  trial  are  plotted  on 
Figs.  14  to  21  .  The  z  coordinate  is  plotted  on  a  logarithmic 
scale  while  the  velocity  u  is  plotted  on  a  linear  scale. 

The  over-water  profiles  were  drawn  as  straight  lines. 

It  can  be  seen  that  not  all  of  the  data  satisfied  the  logarithmic 
"law".  Trials  1,  3,  4,  6  and  10  show  little  deviation  from 
a  straight  line.  However  the  data  points  in  trials  5,  7,  and  8 
are  concave  upwards.  The  profiles  at  an  average  distance  of 
70  meters  downstream  show  either  a  smooth  curve  or  a  straight 
line  in  each  of  the  trials.  This  would  indicate  that  the  wind 
had  approached  equilibrium  with  the  surface,  at  least  in  the  lower 
levels . 


The  intermediate  profiles,  taken  at  an  average  of 
33  meters  from  the  water's  edge,  show  little  if  any  consistent 
pattern.  In  four  of  the  trials  the  velocity  at  the  0.5  meter 
level  was  less  on  the  intermediate  profile  than  in  the  more 
downwind  profile.  All  of  the  intermediate  profiles  have  at 
least  one  point  which  has  a  higher  velocity  than  the  unmodified 
wind  profile.  There  is  no  immediate  explanation  for  this  strange 
behavior.  Anemometers  were  interchanged  both  between  trials  and 
during  trials  in  order  to  eliminate  systematic  error. 

One  feature,  consistent  among  the  intermediate  profiles, 
was  that  the  lower  three  points  on  seven  of  the  eight  trials 
were  nearly  collinear.  It  could  be  said  that  the  layer  from  the 
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surface  to  two  meters  had  adjusted  itself  to  the  new  surface 
conditions.  The  slope  of  each  of  the  intermediate  profiles  is 
seen  to  lie  between  the  slopes  of  the  unmodified  and  the  down¬ 
wind  profiles.  This  indicates  that  the  friction  velocity 
decreased  downwind  from  the  point  of  roughness  change  after 
reaching  a  maximum  value.  The  2  intercept  and  thus  the  roughness 
length  were  different  for  both  of  the  downwind  profiles  in 
each  trial,  except  for  trials  7  and  8.  This  difference  was 
expected  due  to  the  variation  in  the  types  of  grass  encountered 
downwind.  Therefore  one  of  the  conditions  required  by  Townsend’s 
theory  to  allow  self-preserving  flow  has  been  violated. 
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The  modified  Townsend  profile  was  calculated  for  each  of 

the  mean  profiles  using  the  following  method.  In  each  of  the 

four  cases  the  mean  profile  was  first  plotted.  The  values  of 
* 

V  Z1  an<^  Uo  as  obtained  from  the  graphs  were  used  as  the 
parameters  in  the  Townsend  equations.  The  theoretical  profiles 
were  plotted  as  broken  lines  on  Figs.  5,  6,  7,  and  8. 

It  can  be  seen  that  the  Townsend  theory  understimated 
the  degree  of  modification  of  the  profile  in  each  case.  Note 
that  the  slope  of  the  theoretical  lines  runs  parallel  to  that 
of  the  actual  lines  and  also  that  the  roughness  lengths  correspond 
as  would  be  expected. 

In  an  attempt  to  explain  the  discrepancy  between  the 
experimental  and  theoretical  profiles  the  displacement  length  D 
was  introduced.  D  can  be  thought  of  as  a  new  datum  level  below 
which  the  wind  velocity  is  zero.  For  a  first  estimate  D  was 
taken  to  be  25  cm  for  tall  grass.  (Calder  ;  1949).  Then 
each  of  the  profiles  was  also  plotted  for  values  of  D  =  30,  20, 

15  and  10  cm.  For  the  larger  values  of  D  the  profiles 
became  concave  downwards,  while  for  D  =  10  cm  there  was  little 
difference  from  the  case  D  =  0.  Using  the  profiles  for  D  =  15 
the  roughness  lengths  and  friction  velocities  were  measured  and 
these  constants  were  used  to  calculate  the  modified  Townsend 
profiles.  (Figs.  9,10,11,  and  12). 


It  can  be  seen  that  there  was  improved  agreement  between 
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the  theoretical  and  experimental  curves  in  the  neutral  case 
(Fig. 10).  However,  there  was  worse  agreement  in  the  overall 
mean  (Fig.  9)  and  little  difference  in  the  remaining  two  cases 
(Figs. 11  and  12) . 

These  results  serve  to  emphasize  the  importance  of 
using  as  near  neutral  stability  as  possible  in  future 
investigations . 
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CHAPTER  III 


In  Chapter  I  we  looked  at  the  models  presently  available 
which  explain  the  effects  of  a  change  in  surface  roughness.  In 
Chapter  II  we  discussed  and  compared  our  experimental  results  with 
theory.  This  chapter  will  discuss  the  shortcomings  of  the  experiment 
as  conducted  and  offer  suggestions  for  further  investigations. 

In  all  the  theories  that  have  been  published  to  date 
certain  assumptions  have  been  made  which  must  be  recognized  when 
designing  an  experiment.  The  most  important  of  these  assumptions 
is  that  the  atmosphere  is  in  neutral  equilibrium.  Perhaps  in  the 
future  it  will  be  shown  that  to  a  first  approximation  stability 
can  be  ignored.  However,  at  the  present  time  the  majority  of  trials 
should  be  undertaken  in  neutral  conditions.  The  second  assumption 
which  has  been  made  is  that  the  wind  profile  is  fully  established. 

In  other  words  the  shear  stress  is  assumed  constant  with  height. 

Thus  the  second  requirement  is  that  there  is  a  very  long  upwind 
fetch  of  the  order  of  a  kilometer  or  more  (Blom  and  Wartena  1969) . 

A  third  condition  which  must  be  met  is  that  a  single  abrupt  change 
in  roughness  exists  and  the  surface  must  be  of  uniform  composition. 

A  fourth  condition  is  that  the  surface  must  be  nearly  level  so 
as  to  eliminate  systematic  vertical  velocities. 
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The  location  of  the  experiment  is  the  critical  point. 
Several  possibilities  present  themselves.  A  large  body  of 
water  can  be  used  as  the  first  surface  with  a  uniform  shoreline 
as  the  second.  The  advantage  of  this  layout  is  that  an  extremely 
large  discontinuity  can  exist,  the  roughness  length  of  calm 
water  being  of  the  order  0.01  cm  and  that  of  grass  of  the 
order  5  cm.  However,  there  are  a  number  of  disadvantages.  The 
first  is  the  difficulty  of  locating  a  body  of  water  which 
fulfills  the  following  requirements.  The  shore  must  slope  away 
from  the  water  line  very  gently.  The  shore  must  be  of  uniform 
composition  away  from  the  water  for  a  distance  of  100  meters  at 
least.  The  fetch  across  the  water  should  be  of  the  order  of  a 
kilometer.  The  second  disadvantage  is  the  difficulty  encountered 
in  making  measurements  of  wind  and  temperature  over  the  water. 

The  shortcomings  mentioned  above  are  easily  circumvented 
by  using  a  land-to-land  change  of  surface  roughness.  At  Suf field 
the  obvious  choice  for  the  upwind  uniform  surface  is  the  natural 
prairie  grass  which  grows  uniformly  for  tens  of  miles.  Several 
possibilities  present  themselves  for  the  choice  of  the  second 
surface.  The  first  that  comes  to  mind  is  to  use  the  runways  at 
the  airport. 

The  advantages  of  using  runways  are  manifold.  First 
they  are  convenient  and  accessible.  Electric  power  is  nearby  for 
the  instruments  if  needed.  The  orientation  of  the  runways  is  in 
the  direction  of  the  prevailing  winds.  The  airport  receives 
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little  traffic  so  that  there  would  be  little  inconvenience.  In 
the  event  of  rain  the  experiment  won't  be  bogged  down  in  mud. 
There  is  no  danger  of  the  surface  being  trampled  by  people  or 
vehicles  travelling  on  it  as  in  the  case  of  deep  grass.  The 
major  disadvantage  of  the  use  of  a  runway  is  that  the  wind  must 
blow  exactly  parallel  to  it  in  order  to  avoid  lateral  boundary 
effects . 


The  second  thought  that  comes  to  mind  is  to  artificially 
modify  the  prairie  grass  in  order  to  change  its  roughness  length. 
Since  the  grass  -is  naturally  short  it  is  doubtfull  the  roughness 
length  could  be  modified  appreciably  by  mowing.  It  would 
probably  be  more  effective  to  destroy  the  grass  over  a  large  area 
by  either  fire,  herbicide,  or  ploughing.  The  resulting  bare 
earth  would  probably  have  a  roughness  length  appreciably  different 
from  that  of  grass  so  that  the  wind  would  be  significantly 
modified  by  it.  If  such  a  scheme  were  attempted  the  best 
experimental  area  shape  would  be  either  circular  or  octagonal  of 
radius  75  meters.  One  major  advantage  to  such  a  scheme  is  that 
trials  could  be  carried  out  with  winds  from  any  direction. 

A  third  method  of  carrying  out  the  experiment  is  to  use 
two  adjacent  cultivated  fields.  However  such  an  area  is  unlikely 
to  be  found  in  the  vicinity  of  Suffield  where  very  little 
agriculture  is  carried  on. 


My  recommendation  for  a  site  therefore  is  to  modify  the 
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prairie  grass  in  a  flat  region  near  the  station.  The  quantities 
which  are  to  be  determined  are,  primarily,  surface  shear  stress 
and  the  wind  profile  and,  secondarily,  the  temperature  profile. 

In  order  to  make  a  detailed  study  of  the  wind  profiles 
it  is  necessary  to  use  a  large  number  of  levels  and  several  down¬ 
wind  positions.  The  minimum  number  of  levels  I  would  recommend 
is  six,  ranging  from  1/4  meter  to  8  meters  for  the  first 
four  positions;  say  at  10,  20,  40  and  80  meters  from  the 
change  in  roughness  and  from  0.5  to  16  meters  for  downstream 
distances  beyond  80  meters. 

There  will  be  certain  technical  difficulties  associated 
with  such  a  number  of  towers,  so  widely  spaced.  The  first  is  in 
procuring  the  required  materiel;  four  8-meter  towers,  two 
16-meter  towers,  thirty-six  anemometers  and  counters.  Also  there 
would  be  the  sheer  physical  problem  of  aligning  the  towers  with  the 
direction  of  the  wind.  It  was  our  experience  that  two  men 
required  at  least  two  hours  to  set  up  3  towers  and  to  test  the 
anemometers.  Another  major  problem  which  needs  to  be  overcome  is 
that  of  recording  all  the  data  simultaneously.  To  connect  all 
anemometers  by  cable  to  a  central  point  would  require  close  to 
500  meters  of  7— conductor  cable;  (or  13-  conductor  cable  if 
temperature  measurements  are  also  made) . 

These  problems  are  such  that  it  would  be  mandatory  to  set 
the  towers  up  permanently  and  then  wait  for  the  wind  to  blow  in  the 
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chosen  direction.  In  order  to  overcome  these  difficulties  without 
losing  resolution,  a  single  highly  portable  tower  could  be 
employed  like  a  probe  in  a  wind  tunnel.  It  could  be  mounted  on 
the  rear  of  a  half  ton  truck,  together  with  the  counters,  and  be 
moved  from  point  to  point  downwind. 

It  is  necessary  to  assume  a  steady  state  wind  when 
using  a  moveable  tower.  This  condition  is  satisfied  to  a  high 
degree  whenever  the  geostrophic  wind  is  relatively  strong. 

The  -experimental  technique  would  be  as  follows.  One 
tower  would  be  set  up  ahead  of  the  roughness  change  and  to  one 
side.  -  Ideally  it  would  be  a  16  meter  tower  but  an  8  meter 
tower  would  suffice.  One  man  would  be  required  to  record  the 
readings  on  counters  every  15  minutes.  The  moveable  probe  would 
be  positioned  first  at  10  meters  away  from  the  line  of  change 
and  a  ten  minute  mean  profile  obtained.  These  profiles  might  be 
taken  at  10,  20,  40,  80  or  any  desired  interval  depending  upon 
what  type  of  cross-section  is  desired. 

Since  the  equipment  is  highly  portable  it  is  necessary  to 
be  able  to  -measure  distances  accurately  and  quickly.  A  suggestion 
would  be  -to  count  the  number  of  revolutions  of  the  wheels  of  the 
truck,  and  knowing  their  circumference  to  calculate  the  distance 
travelled.  The  most  practical  way  to  count  wheel  revolutions  is 
to  use  the  truck’s  speedometer  cable.  A  cam  on  the  cable  can  be 
used  to  actuate  a  microswitch.  A  standard  electric  counter  can 
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then  be  used.  Normally  the  cable  makes  1000  revolutions  per 
mile  or  1  revolution  per  5.28  feet. 

A  tower  or  mast  should  be  attached  to  the  rear  of  a 
pickup  truck  in  such  a  way  as  to  allow  it  to  be  lowered  for 
travelling  and  for  making  adjustments  to  the  instruments.  A 
practical  method  of  building  the  tower  on  a  truck  is  shown  in 
Figure  13  .  The  mast  should  be  constructed  of  two-inch  steel 
tubing.  Anemometer  supports  would  be  welded  in  place  so  as  to 
place  the  cups  at-  8.0,  4,0  and  2.0  meters. 

-A  frame  constructed  of  angle  iron  is  bolted  to  the  rear 
of  the  truck  and  places  the  mast  pivot  slightly  above  the  level 
of  the  cab.  Two  guy  wires  are  attached  to  the  sides  of  the  body 
of  the  truck.  To  eliminate  the  need  for  a  third  rope  to  the 
ground,  an  aluminum  pole  of  2  inch  diameter  is  used.  It  is 
bolted  to  a  bracket  welded  to  the  mast  at  roughly  the  six  meter 
level  and  is  carried  in  a  socket  on  the  rear  edge  of  the  cab. 

The  anemometers  for  the  2„0,  1.0,  0.5  and  0.25  levels 
are  mounted  on  a  separate  stand  constructed  of  water  pipe.  This 
allows  the  stand  to  be  located  several  meters  to  one  side  of  the 
truck.  The  lower  level  instruments  therefore  will  be  outside 
the  area  where  -the  wind  would  be  disturbed  by  the  presence  of  the 
truck. 


The  following  sketches  will  indicate  a  possible  method  of 


constructing  the  apparatus:  (Fig.  13) 
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SUMMARY 


The  object  of  this  project  was  to  find  experimental 
confirmation  of  one  of  the  theories  describing  the  changes  in 
wind  profile  and  associated  parameters  to  a  change  in  surface 
roughness  under  conditions  of  neutral  stability. 

For  a  number  of  reasons  which  have  been  discussed 
previously,  the  desired  quantitative  results  could  not  be 
achieved . 

However,  a  few  qualitative  conclusions  could  be  drawn 
from  visual  examination  of  the  wind  profiles.  First  the  lowest 
layer  of  air  (  O  -  1.0  meter)  appears  to  come  to  equilibrium 
with  a  new  surface  very  quickly  and  is  closely  described  by  the 
logarithmic  profile,  thus  supporting  one  of  the  assumptions  of 
Townsend.  Second  the  friction  velocity  changes  abruptly  at  or 
close  to  the  point  of  change  of  roughness.  It  appears  to  return 
to  a  new  equilibrium  value  far  downstream.  Third  for  one  case 
which  was  assumed  to  be  under  conditions  of  near  neutral  stability 
the  Townsend  profile  fits  the  data  closely.  However,  the  lack  of 
accurate  stability  measurements  prevents  any  firm  conclusions  to 


be  drawn. 
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APPENDIX  A 


(Elliott  1958) 


THE  ELLIOTT  PROFILE 

The  problem  Elliott  undertook  to  solve  was  to  compute 
the  height  of  the  internal  boundary  layer  as  a  function  of 
downwind  distance  x.  First  he  considered  a  region  bounded  by 
the  top  of  the  boundary  layer,  the  ground  and  two  vertical  lines 
a  distance  x  apart.  The  loss  of  momentum  out  of  the  region  by 
convection  according  to  the  Von  Kdrm&n  integral  theorem  is  given 
by 


2,  d 
PU  dz  -  U,  -r~ 
h  dx 


pudz}6x. 


(Al) 


Assuming  conservation  of  momentum  this  expression  must  equal  the 
net  gain  of  momentum  due  to  vertical  flux  into  the  region.  The 
vertical  momentum  flux  is  equal  to  the  difference  in  shearing 
stress  between  the  top  and  the  bottom  of  the  layer.  Thus 

-y-  J  pu^dz  -  u,  I  Pudz  =  x  -  x  .  (A2) 

dx  J  h  dx  J  ho 

z  z 

o  o 

This  equation  can  now  be  used  to  determine  h  if  the  variation 
of  wind  speed  with  height  is  known. 
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'"tt 

Elliott  assumed  a  neutral  atmosphere  with  u*  =  —  In  z/z' 

k  o 

in  the  undisturbed  region  upwind  of  and  above  the  disturbed  area. 

The  wind  profile  in  the  internal  boundary  layer  was  assumed  to  be 
u* 

u  =  In  z/zo.  When  these  profiles  are  used  it  is  implied  that 
uj.  is  constant,  because  of  equilibrium  conditions  upstream,  and 
that  u^  is  a  function  of  x  but  not  of  z.  Therefore,  at  the 
top  of  the  boundary  layer  there  will  be  a  discontinuity  in  the 
stress  and  in  the  velocity  gradient.  With  these  assumptions 
equation  (A2)  was  written  as: 


_d_ 

dx 


h 

/ 


z 

o 


2, 

u  dz 


d 

Uh  dx 


h 

/ 


z 

o 


udz 


(A3) 


The  full  solution  of  (A3)  was  derived  by  Elliott  in  his  Ph.D. 
dissertation  in  1958.  He  obtained  the  following  differential 
equation: 


{ln  zi  - 3  -  ^  + 1 rz 


l  ^Z  In  Z  2 

7  (1~  z“)}“dT  =  (1  +  i^TSz7)r_ 

11  1  o 


where  Z=h/z  ,  m=z/z’  . 
1  o  oo 


Integrating  (A4)  results  in 


00  (In  z  )n 

2kX  =  Z  In  Z  -  4Z  -  In  Z  +  4  l  - = — 

111  1  L.  n*n: 

n=l 


ln(ln  cZ^) 


+  In  c{ 


c 


(In  Zj^-3-20) 
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,  °°  (In  cZ1  )n 

+  C  l  <ln  Z1  -  3» 

n=l 


i  ~  u»  v 

+  4 (In  In  cZ  )(1  -  ■=■)  +  4-{  T  - — 

J-  c  n  •  n . 

n=l 


n 


n 


00  (In  cZ  ) 

-  ~  /  - i - }  +  constant. 

c  n  •  n . 

n=l 


The  constant  c  was  determined  by  the  condition  that 
h  =  Zq  at  x  =  0  which  may  not  be  true  in  practice.  From 
dimensional  arguments  Elliott  derived  an  approximation  to  his 
solution: 


(7b  (In  -7- 
z  z 

o  o 


1/5 
)  a  X 


1/5 


(A5) 


APPENDIX  B 


(Panofsky  and  Townsend  1964) 


PANOFSKY  AND  TOWNSEND  PROFILE 

This  theory  is  an  extension  of  the  Elliott  theory  and 
differs  in  only  one  detail.  Whereas  Elliott  assumed  a  logarithmic 
profile  in  the  disturbed  region,  Panofsky  and  Townsend  used  a 
log-linear  profile. 

* 

u.. 

u  =  -y  [  (1-s)  In  —  +  s  |]  (Bl) 

o 

*  *  * 
where  s  =  (u,  -  u  )/u. . 

101  * 

The  profile  follows  from  the  mixing  length  relation  -r—  =  — 

o  z  kz 

with 


*  *  z ,  * 

U  =  U  [  (1-s)  +  S  -r]  =  u  + 
1  do 


*  *  7. 

<vuo>! 


from  which  it  can  be  seen  that  the  friction  velocity  is  a  linear 
function  of  height  with  the  value  u*  =  u*  at  z/d  =  0  and 
u*  =  u*  at  z/d  =  1.  Choose  a  streamline  at  a  level  z^  well 
upstream  from  the  discontinuity.  Then  the  volume  of  air  per  unit 
time  passing  a  vertical  plane  is 


/ 


udz 


ZjJln^/ z^)  -  1]  . 


(B2) 


The  height  of  the  same  streamline  at  the  boundary  is  d  and  the 
volume  of  air  is 
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d 

J  udz 
z 

o 


[(l-s)(ln(d/zo)  -1)  +  s/2]d 


(B3) 


Assuming  steady  and  incompressible  flow  these  expressions  are  equal. 
Since  the  velocity  must  be  continuous  across  the  interface 


(l-s)ln(d/z  )  +  s  =  In  z  /z'  .  (B4) 

o  1  o 


Eliminating  ^ 

1  -  (ln(d/z  ) (1- 

s (ln(d/z  )-l)  =  In  z'/z  -  ln[ - 5 - - ] 

°  °°  1  -  (ln(d/ zQ))-1 

(E5) 


in  which  the  last  term  is  negligible  ip  most  cases.  The  approximation 


In  z'/z 
_ o  o 

ln(d/z  )-l 
o 


(B6) 


is  sufficiently  accurate  and  implies  that  there  is  no  displacement 
of  the  streamlines  caused  by  the  discontinuity  in  roughness. 


Conservation  of  momentum  in  the  whole  flow  is  imposed  by 
the  integral  constraint: 


t  ,  3u  3u.  ,  *2 

J  (U  3?  +  W  H)dZ  "  U1 


*2  *2 

u  =  u,  s(2-s) 
o  1 


(B7) 
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appendix  c 


(Taylor  1962) 


THE  TAYLOR  THEORY 

Taylor  considered  the  situation  as  sketched  below: 


A  fully  extablished  profile  in  Region  1  encounters  an 
increased  roughness  at  the  origin,  point  0.  In  Region  2,  a 
new  profile  with  constant  shearing  stress  is  established.  The 
region  OBDC  is  the  region  of  transition  with  0C  and  BD  as 
lines  of  constant  stress.  The  line  PQ  represents  a  typical 
streamline.  The  length  between  the  origin  and  a  point  midway  in 
the  transition  region  is  called  n  .  Taylor  derived  an  expression 
for  n  as  a  function  of  z. 


Taylor's  starting  point  was  the  equation  of  motion: 


47. 


Defining  a  stream  function  <j>  ,  with 

Tz“  =  u’  IT'  “  "  w  » 

he  rewrote  the  equation  of  motion 


(C2) 


^1  3_u  _  _^1  3u  _  _8_  2  _  _1  g>£ 

9 z  9x  9x  9z  9 z  U*  p  3x  * 


(C3) 


Integration  over  the  transition  region  gave: 


f  (u"-u’)dp1  =  -  f  u^dx  “  P  1  £  f>  ^ (C4) 
o 


=  -  X  -  Y  . 


Single  and  double  primes  represent  quantities  in  Regions  1  and  2 
respectively.  <)>  is  the  value  of  the  stream  function  along  CD. 

The  velocity  difference  (u"  -  u')  is  measured  along  a  streamline. 
The  line  integrals  are  taken  in  the  direction  OBDC.  X  was 
shown  to  be: 


n(u^2  -  u”2) 

and  Y  was  shown  to  be  small  compared  with  X.  The  stream  function 


for  wind  in  neutral  conditions  was  shown  to  be: 


■  ^  ~ 

? 
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<P  =  — {  ln(z/ z^)  -  1} 


(C5) 


for  z/z  >>  1  and 


o 


(C6) 


o 


It  is  now  possible  to  calculate  n  given  u'  ,  u"  ,  z'  and  z" 

*  *  o  o 

using  equations  (C4) ,  (C5)  and  (C6)  without  performing  any 
integrations. 

By  means  of  dimensional  analysis  and  by  an  examination  of 
the  available  data,  Taylor  came  to  the  conclusion  that  the  ratio 
of  the  new  friction  velocity  to  the  old  was  equal  to  the  ratio 
of  the  roughness  lengths  raised  to  the  0.09  power.  From  wind 
tunnel  experiments  he  concluded  that  the  friction  velocity  was 
constant  downwind  from  the  change  of  roughness. 
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appendix  d 


TRIAL  1 

18  July  68 

Wind  Direction  205° 

Orientation  of  Towers  185° 

Cloud  Cover  none 

Temperature  79°F. 

Location  D 

Fetch  350  m. 

Water  Temperature  80° F. 


x  =  3 

33 

63 

1 

8.0 

3.10 

3.10 

3.04 

4.0 

2.90 

2.97 

2.86 

2.0 

2.71 

2.67 

2.36 

1.0 

2.55 

2.17 

1.93 

0.5 

2.09 

1.51 

1.51 

(x  and  z  measured  in  meters;  other  numbers  are  velocity  in  m/sec.) 


' 
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TRIAL  2 

31  July 

Wind  Direction 

Orientation  of  Towers 

Cloud  Cover 

Temperature 

Location 

Fetch 

Water  Temperature 


215° 
220° 
none 
70°F . 
D 

350m. 
70°F . 


x  =  0 

6 

16 

8.0 

5.85 

6.26 

4.0 

5.83 

5.75 

2.0 

5.34 

5.60 

5.44 

1.0 

5.26 

4.31 

0.5 

4.61 

3.10 

2.35 

TRIAL  3 

12  August 
Wind  Direction 
Orientation  of  Towers 
Cloud  Cover 
Temperature 
Location 


330° 

310° 

100% 

63°F. 

A 


Fetch  280m. 

Water  Temperature  66°F. 


51 


x  =  3 

19 

35 

8.0 

7.52 

7.39 

7.35 

4.0 

7.07 

7.15 

6.86 

2.0 

6.28 

6 . 56 

5.93 

1.0 

5.52 

5.37 

4.60 

0.5 

4.93 

3.73 

3.35 

TRIAL  4 

13  August 


Wind  Direction  100° 

Orientation  of  Towers  90° 


Cloud  Cover 
Temperature 
Location 
Fetch 

Water  Temperature 


scattered  Cu 
67°F. 

B 

300m . 

68°F. 


z\ 

x  =  3 

33 

78 

8.0 

5.73 

5.99 

5.77 

4.0 

5.35 

5.64 

4.86 

2.0 

5.00 

4.89 

4.35 

1.0 

4.43 

3.86 

3.70 

0.5 

3.95 

2.75 

2.92 

52 


TRIAL  5 

14  August 

Wind  Direction 

Orientation  of  Towers 

Cloud  Cover 

Temperature 

Location 

Fetch 

Water  Temperature 


90° 

90° 

Overcast,  occasional  drizzle 
54°F. 

B 

300m. 

58°F . 


z\ 

x  =  3 

33 

78 

8.0 

5.69 

5.42 

5.60 

4.0 

5.26 

5.57 

5.12 

2.0 

4.91 

4.72 

4.20 

1.0 

4.43 

2.74 

3.41 

0.5 

3.75 

2.53 

2.83 

TRIAL  6 


16  August 

Wind  Direction 

Orientation  of  Towers 

Cloud  Cover 

Temperature 

Location 

Fetch 

Water  Temperature 


90° 

90° 

overcast  thin  cirrus 

64°F . 

C 

300m . 

62°F . 


f 
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z  \ 

x  =  3 

33 

63 

8.0 

3.10 

3.21 

3.04 

4.0 

2.97 

2.92 

2.79 

2.0 

2.38 

2.47 

2.39 

1.0 

2.36 

1.83 

1.85 

0.5 

2.22 

1.00  ! 

1.24 

TRIAL  7 

18  August 

Wind  Direction 

Orientation  of  Towers 

Cloud  Cover 

Temperature 

Location 

Fetch 

Water  Temperature 


330° 

330° 

broken  Cu 
64°F . 

A 

280m. 

67°F . 


l\ 

x  =  3 

13 

23 

8.0 

l 

7.33 

7.79 

7.34 

4.0 

7.10 

6.99 

6.79 

2.0 

6.54 

6.71 

5.83 

1.0 

5.97 

4.76 

4.10 

0.5 

5.15 

3.05 

2.38 
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TRIAL  8 

20  August 

Wind  Direction 

Orientation  of  Towers 

Cloud  Cover 

Temperature 

Location 

Fetch 

Water  Temperature 


z  \ 

x  =  3 

13 

8.0 

6.31 

6.58 

4.0 

6.13 

6.03 

2.0 

5.82 

5.83 

1.0 

5.42 

4.32 

0.5 

4.27 

2.65 

TRIAL  9 

26  August 

Wind  Direction 

Orientation  of  Towers 

Cloud  Cover 

Temperature 

Location 

Fetch 


90° 

90° 

overcast,  light  rain 
63°F. 

C 

300m. 

65°F. 

23 

6.23 

5.66 

4.90 

3.54 

2.10 


210° 
190° 
8/10  As 

E 

350m. 


Water  Temperature 


’ 
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2  \ 

x  =  4 

34 

64 

8.0 

10.09 

10.39 

8.58 

4.0 

9.63 

9.56 

8.95 

2.0 

8.78 

9.03 

7.44 

1.0 

8.30 

6.42 

5.49 

0.5 

7.54 

4.48 

3.53 

0 


5  6  . 


U  (m/s) 


F,3  IS 


~J~<r  i  a  I  ^ 


si. 


U.  (  nr\  /  s) 

Fl3  H  “Trial  5 


—  — 
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'  V£  -  ^ 


X 


\ 
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F13  n  Tr<«/  7 


S  9. 


Z  (m) 


u.  m  /  S) 
Trlal  9 
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60. 


OVERALL  MEAN  PROFILE 
(using  data  from  Trials  1,  4,  5,  6,  9) 


z\ 

x  =  3.0 

33.0 

69 

8.0 

5.55 

5.62 

5.21 

4.0 

5.22 

5.33 

4.92 

2.0 

4.76 

4.76 

4.15 

1.0 

4.42 

3.41 

3.28 

0.5 

3.91 

2.45 

2.40 

MEAN  OF  "NEUTRAL"  TRIALS 

(5,9) 

z  \ 

x  =  3.5 

33.5 

71 

8.0 

8.42 

8.44 

7.57 

4.0 

7.95 

8.08 

7.51 

2.0 

7.31 

7.34 

6.21 

1.0 

6.80 

4.89 

4.75 

0.5 

6.02 

3.74 

3.39 

MEAN  OF  UNSTABLE 

TRIALS  (1,4,6) 

2  \ 

x  =  3.0 

33.0 

68.0 

8.0 

4.25 

4.38 

4.22 

4.0 

3.99 

4.10 

3.74 

2.0 

3.59 

3.57 

3.24 

1.0 

3.32 

2.80 

2,66 

0.5 

2.94 

1.87 

2.02 

61. 


MEM  OF  TRIALS  7  MD  8 


z  \ 

x  =  3 

13 

23 

8.0 

6.81 

7.18 

6.79 

4.0 

6.61 

6.50 

6.20 

2.0 

6.19 

6 . 26 

5.35 

1.0 

5.70 

4.54 

3.82 

0.5 

4.71 

2.85 

2.25 

